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Abstract 

We have derived the interaction between Jahn- Teller centers due to optical and acoustic phonons. 
Without Coulomb repulsion the model leads to a global phase separation at low temperature. At 
long distances Coulomb repulsion always dominates the short range attraction leading to phase 
separation on a short scale. On the basis of Monte-Carlo simulation of microscopic lattice-gas 
model we have formulated an effective phenomenological model. In the limit of low density of 
polarons, phase separation takes place in the form of charged bubbles. In the limit of high density, 
charge segregation occurs on domain walls. 
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There is substantial evidence that the ground state in many of oxides is inhomogeneousQ]. 
In cuprates, for example neutron-scattering experiments surest that phase segregation takes 
place in the form of stripes or short segments of stripes^, O]. There is some controversy 
whether this phase segregation is associated with magnetic interactions. On the other hand 
it is also generally accepted that the charge density in cuprates is not homogeneous. 

The idea of charge segregation in cuprates appeared just after discovery of High-Tc 
materials 0,0,1^. In most cases long-range Coulomb repulsion was not considered. Recently 
it was suggested that interplay of the short range lattice attraction and the long-range 
Coulomb repulsion between charge carriers could lead to the formation of short metallic M, ^1 
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or insulating [8|,|9|| stripes of polarons. If the attractive potential is isotropic, charged bubbles 
have a spherical shape. Kugel and Khomskii suggested recently that the anisotropic 
attraction forces caused by Jahn- Teller centers could lead to the phase segregation in the 
form of stripes. The long-range anisotropic attraction forces appear as the solution of the full 
set of elasticity equations (see ref.jlll). Alternative approach to take into account elasticity 
potentials was proposed in ref. and is based on the proper consideration of compatibility 
constraint caused by absence of a dislocation in the solid. 

Rece,.>. we fonnu.a.ed .ode.Q whe.e we suggested ,„te..ac«o„ of a .wo-fo.d 
degenerate electronic state with ti phonon modes at a finite wave-vector can lead to a 
ocal nonsymmetric deformation and short-length scale charge segregation. It was suggested 
1^ that this phase separation may lead to the insulator to superconductor phase transition 
governed by percolation. In this paper we reduce the proposed model to the lattice gas model 
and show that the model indeed displays phase separation, which may occur in the form of 
stripes or clusters depending on anisotropy of the short range attraction between localized 
carriers. We also generalize the model taking into account interaction of the Jahn- Teller 
centers via elasticity induced field. 

We first derive a simplified lattice gas model. We show that the model without Coulomb 
repulsion displays a first order phase transition at a constant chemical potential. When 
the number of particles is fixed, the system is unstable with respect to the global phase 
separation below a certain critical temperature. In the presence of the Coulomb repulsion, 
global phase separation becomes unfavorable due to a large contribution to the energy from 
long range Coulomb interaction. The system shows mesoscopic phase separation where the 
size of charged regions is determined by the competition between the energy gain due to 
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ordering and energy cost due to breaking of the local charge neutrality. Since the short 
range attraction is anisotropic the phase separation may be in the form of short segments 

or /and stripes. ^ ^ 

Let us begin by considering a simplified version of the JT model Hamiltonian taking 
only the deformation of the Big symmetry: 

i^JT = ^7E^3,l/(n)(6l+„ + 6l+n), (1) 
n,l 

here the Pauli matrix a^^i describes two components of the electronic doublet, and /(n) = 
(n^ — n'^)fo{n) where /o(^) is a symmetric function describing the range of the interaction. 
For simplicity we omit the spin index in the sum. The model could be easily reduced to a 
lattice gas model. Let us introduce the classical variable ^i=<hf + h-^> / \pl and minimize 
the energy as a function of $i in the presence of the harmonic term c^X^i $^/2. We obtain 
the deformation, which corresponds to the minimum of energy, 

$S°) = -v^(7/a;E^3,i+n/(n). (2) 

n 

Substituting ^P"* into the Hamiltonian (1) and taking into account that the carriers are 
charged we arrive at the lattice gas model. We use a pseudospin operator S* = 1 to describe 
the occupancies of the two electronic levels niand Here = \ corresponds to the state 
with rii = 1 , n2 = 0, SI = —1 to rii = 0, n2 = 1 and = to rii = n2 = 0. Simultaneous 
occupancy of both levels is excluded due to a high onsite Coulomb repulsion energy. The 
Hamiltonian in terms of the pseudospin operator is given by 

Hlg = - J)^f^/ + K(i - j)QiQj), (3) 

where Qi = {S?Y. Vc{n) = e^/eoa(n^ + riyY^'^ is the Coulomb potential, e is the charge 
of electron, eo is the static dielectric constant and a is the effective unit cell period. The 
anisotropic short range attraction potential is given by Vi{n) = g'^/ujJ2m /(iii)/(ii+ni). The 
attraction in this model is generated by the interaction of electrons with optical phonons. 
The radius of the attraction force is determined by the radius of the electron-phonon inter- 
action and the dispersion of the optical phonons 8| . 

A similar model can be formulated in the limit of continuous media. In this case the 
deformation is characterized by the components of the strain tensor. For the two dimensional 
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case we can define 3 components of tlie strain tensor: ei = u^x + Uyy, e = Uxx — Uyy and 
62 = Uxy transforming as the Aig, Big and B2g representations of the D^h group respectively. 
These components of the tensor are coupled linearly with the two-fold degenerate electronic 
state which transforms as the Eg or representation of the point group. Similarly to 
the case of previously considered interaction with optical phonons we keep the interaction 
with the deformation e of the Big symmetry only. The Hamiltonian without the Coulomb 
repulsion term has the form: 

H = gT. S-ei + \ {Aiel, + A.ej + ^36^;) , (4) 

i 

where Aj are corresponding components of the elastic modulus tensor. The components of 



the strain tensor are not independent 



12l lisj l and satisfy to the compatibility condition. 
The compatibility condition leads to a long range anisotropic interaction between polarons. 
The Hamiltonian in the reciprocal space has the form: 

^ = ^E^kek + (A2 + Ai[/(k))|. (5) 
The Fourrier transform of the potential is given by: 

^^^^ = FTWfe^- 

By minimizing the energy with respect to and taking into account the long-range 
Coulomb repulsion we again derive Eq.(3). The anisotropic interaction potential VJ(n) = 

2 

— I^kGxp (^k ■ n) 2(A2+Ait/(k)) determined by the interaction with the classical deformation 
and is long-range. It decays as at large distances in 2D. Since the attraction forces de- 
cay faster then the Coulomb repulsion forces at large distances the attraction can overcome 
the Coulomb repulsion at short distances, leading to a mesoscopic phase separation. 

Irrespective of whether the resulting interaction between polarons is generated by acoustic 
or optical phonons the main physical picture remains the same. In both cases there is 
an anisotropic attraction between polarons on short distances. The interaction could be 
either ferromagnetic or antiferromagnetic in terms of the pseudospin operators depending 
on mutual orientation of the orbitals. Without loosing generality we assume that V{y\) is 
nonzero only for the nearest neighbors. 

Our aim is to study the model (Eq.(3)) at constant average density, 

^=^EQi' (7) 



where A^is the total number of sites. However, to clarify the physical picture it is more 
appropriate to perform calculations with a fixed chemical potential first by adding the term 
—f^J2iQi to the Hamiltonian (3). 

Models such as (3), but in the absence of the long-range forces, were studied many years 



ago on the basis of the molecular-field approximation in the Bragg- Williams formalism 
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17l |. The mean-field equations for the particle density n and the pseudospin magnetization 
M= ^Ei S-^ have the formQ: 

^ ^ 2 sinh {2zViM/knT) 

exp {-fi/k^T) + 2 cosh {2zViM/kBT) ^ ' 

^ 2 cosh {2zViM/kBT) 

^ exp {- 12/ keT) + 2 cosh {2zViM/kBT) ^ ^ 

here z = 4 is the number of the nearest neighbours for a square lattice in 2D and ks is 
the Boltzman constant. A phase transition to an ordered state with a finite M may be of 
either the first or the second order, depending on the value of fi. For /i > it is always of 
the second order and n — 1 as T ^ 0. For the large negative values fi < —2zVi the phase 
transition is absent and n ^ as T ^ 0. For the physically important case ~2zVi < fi < 
ordering occures as a result of the first order phase transition. Two solutions of Eqs.(8,9) 
with M = and with finite M corresponds to two different minima of the free energy. The 
line of the phase transition is determined by the condition: F{M = 0, fi, T) = F{M, fi, T) 
where M is the solution of Eq. (8) (Fig.l)). When the number of particles is fixed (Eq.9) 
the system is unstable with respect to global phase separation below Tcru- As a result at 
fixed average n two phases with uq = n{M = 0, /i, T) and um = n{M, yU, T) coexist. 

To investigate further the properties of the system, we performed Monte-Carlo (MC) sim- 
ulations of the Hamiltonian Eq.(3). The simulations were performed on a square lattice with 
dimensions up to L x L sites with 10 < L < 100 using a standard Metropolis algorithm [isj] 
in combination with simulated annealing [iQ^. 

First, for comparison with MF theory, the Monte-Carlo simulation of the model Eq.(3) 
in absence of the Coulomb forces shows the reduction of Tcru due to fluctuations in 2D, by 
a factor of ~ 2 (Fig.l). 

Next we include the Coulomb interaction Vc{r). We use open boundary conditions to 
avoid complications due to the long range Coulomb forces and ensure global electroneutrality 
by adding an electrostatic potential due to the uniformly charged background Vjeii (i) to 
(3), such that the Monte Carlo interaction becomes Hmc = H^q + J2iQiVjeii (i), with the 



total background charge being equal in magnitude to the total charge of carriers, eJ^iQu 
but with opposite sign. The dimensionless temperature t = fc^Teoa/e^, and the energy per 
particle are defined as cmc = Hmc^qO^/^^^- The short range potential = V;(i)eoa/e^ 
was taken to be nonzero only for |i| < 2 and was therefore specified for nearest neighbours 
and next nearest neighbours as vi{l,0) and vi{l, 1) respectively. 

The results of the Monte Carlo simulations could be summarized as follows, (i) Coulomb 
repulsion leads to further reduction of the onset temperature of phase separation (see Fig.l). 
(m) Depending on the value of the attractive potential and its anisotropy phase separation 
takes place in the form of bubbles, horizontal and vertical or diagonal stripes (Fig. 2). {in) 
Depending on the anisotropy of the potential stripes (bubbles) are ferromagnetically or 
antiferromagnetically orbitally ordered (Fig.2). (iv) The size of clusters is determined by 
the ratio of short range attraction and long-range Coulomb energy. 

The results of the Monte-Carlo simulation of the model (3) allow general model indepen- 
dent interpretation. Let us consider the classical free energy density corresponding to the 
first order phase transition: 



Here t is dimensionless temperature. Let us assume that the order parameter A is coupled 
to the local charge density p, Fcoupi = a(A^ ~ pY- The total free energy density should 
contain the gradient term Fgrad = C(VA)^ and electrostatic energy F^i = K(j)p as well. The 
electrostatic potential is determined from the Poisson equation eoA0 = 47re(p — po)- Here 
we write po explicitly to take into account global electroneutrality. Substituting the solution 
of the Poisson equation to the electrostatic energy we obtain Fei = K' Y.kyik){p — po)k5 
here V{k) = 1/k. Minimization of Fi + Fcoupi + Fgrad + F^i at fixed t and po gives the spatial 
variation of the order parameter. The results of minimization are presented in Fig. 3. As 
clearly seen from Fig. 3 at low density phase separation takes place in the form of bubbles 
of optimal radius which are ordered. When density increases charge segregation appears in 
the form of charged domain walls. 

Let us assume that a single bubble of the ordered phase with the radius R has appeared. 
As it was discussed in jl3l.]2c| the energy of the bubble is determined as: e = —FnR^ + 
anR + '-fuR^, here F is the energy difference between two minima in the free energy, a is 
the surface energy, and 7 determines the energy due local breaking of the charge neutrality. 




(10) 
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li a < ttF/S^ e has well defined minimum. This minimum corresponds to optimal size of 
the cluster. Residual interactions leads to interbubble interactions and ordering of clusters 
at low temperatures. 

We have demonstrated that anisotropic interaction between Jahn- Teller centers generated 
by optical and/or acoustical phonons leads to the short scale phase separation in the presence 
of the long range Coulomb repulsion. Topology of texturing differs from charged bobbles to 
oriented charged stripes depending on the anisotropy of short range potential. 
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I. FIGURE CAPTIONS 

Fig. 1 The phase diagram of the model in absence of the Coulomb repulsion. The full 
symbols represent the mean field (MF) solution while empty symbols represent the Monte- 
Carlo (MC) simulation result on systems with two different sizes L. 

Fig. 2 Snapshots of clusters ordering at t = 0.04, n = 0.2 and ^^(l, 0) = — 1 as a function 

of Vi{l,l). 

Fig. 3 Snapshots of the order parameter distribution at different average densities ob- 
tained by minimization of the free energy (10). With increasing of density crossover from 
charged bubbles to the charged domain walls is clearly observed. 
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